Transient absorption and reshaping of ultrafast xuv light by laser-dressed helium 
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We present a theoretical study of transient absorption and reshaping of extreme ultraviolet (xuv) 
pulses by helium atoms dressed with a moderately strong infrared (ir) laser field. We formulate the 
atomic response using both the frequency-dependent absorption cross section and a time-frequency 
approach based on the time-dependent dipole induced by the light fields. The latter approach can 
be used in cases when an ultrafast dressing pulse induces transient effects, and/or when the atom 
exchanges energy with multiple frequency components of the xuv field. We first characterize the 
dressed atom response by calculating the frequency-dependent absorption cross section for xuv 
energies between 20 and 24 eV for several dressing wavelengths between 400 and 2000 nm and 
intensities up to 10 12 W/cm 2 . We find that for dressing wavelengths near 1600 nm, there is an 
Autler-Townes splitting of the Is — >• 2p transition that can potentially lead to transparency for 
absorption of xuv light tuned to this transition. We study the effect of this xuv transparency in a 
macroscopic helium gas by incorporating the time- frequency approach into a solution of the coupled 
Maxwell-Schrodinger equations. We find rich temporal reshaping dynamics when a 61 fs xuv pulse 
resonant with the Is — > 2p transition propagates through a helium gas dressed by an 11 fs, 1600 nm 
laser pulse. 

PACS numbers: 32.80.Fb, 32.80.Qk, 32.80.Rm, 42.50.Hz 



I. INTRODUCTION 

The advent of ultrafast xuv and even x-ray light 
sources that can be synchronized to optical or ir laser 
pulses has given rise to several recent studies of the tran- 
sient absorption of such radiation by laser-dressed atoms, 
both experimental p]-[5| and theoretical . For exam- 
ple, many of the experiments done in attosecond physics 
involve the transient absorption of attosecond xuv radi- 
ation by atoms interacting with an ir laser field. This is 
because the strong-field process of high harmonic gen- 
eration (hhg), which is used to produce the attosec- 
ond xuv radiation as either single pulses or trains of 
pulses, results in the xuv field being precisely synchro- 
nized with the driving IR field [H, Q. It is then possi- 
ble to perform experiments using the xuv field and a 
replica of the original IR field with attosecond precision 
jlOj . Glover et al. also showed that it is possible to over- 
lap pulses of synchrotron-produced soft x-ray radiation 
with an 800 nm dressing laser in a study of laser-induced 
transparency in neon |2|. 

In this paper we explore how an ultrafast xuv pulse 
interacts with a simple atom, helium, in the presence of 
a moderately strong ir field which may be either shorter 
or longer in duration than the XUV pulse. We have as our 
goal formulating theoretical methods that can be used to 
calculate the absorption and emission of xuv radiation 
by strongly dressed atoms even when the xuv pulses are 
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on the femtosecond time scale, and may include multiple 
frequencies in a comb. In addition we want to be able to 
study the transient absorption and reshaping of radiation 
as it propagates through a macroscopic amount of gas. 
We will restrict ourselves in this study to cases where the 
IR laser dresses the atom without appreciably exciting it, 
leaving higher IR intensities for a future paper. We will 
also restrict ourselves to xuv wavelengths and intensities 
where single excitations below the first ionization thresh- 
old at 24.6 eV dominate the XUV absorption. Even given 
these restrictions, the IR laser has a substantial impact 
on the Rydberg and continuum states of the atom and, 
in this way, enab les p rofound control over resonant XUV 
absorption @, Ill14l4l]. 

In the calculations we present we will consider the sim- 
plest case, where there are just two radiation fields, one 
which dresses the atom and one which is absorbed and 
possibly reshaped. The fundamental problem of a two- 
color field like this has been studied before in the con- 
text of x-ray absorption by neon 0, 0, [H-[I3 > argon [l8[ , 
and krypton (l9| atoms. Specifically, the examination of 
laser-dressed atoms led to the discovery of electromag- 
netically induced transparency (eit) for x-ravs 0], bet- 
ter characterized as Autler-Townes splitting [12j because 
the transparency is not predominantly caused by destruc- 
tive interference. There have also been several studies of 
helium in the context of the two-color problem we are 
discussing. It was investigated with an optical laser and 
the xuv free electron laser in Hamburg |20rl22l |. and the 
impact of laser-dressing helium on the production of xuv 
radiation via high harmonic generation (hhg) was stud- 
ied theoretically in 23] . 

We begin by characterizing the single atom response 
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in terms of the cross section for absorption of xuv ra- 
diation of frequency u)x- First we calculate the linear, 
frequency-dependent xuv absorption cross section using 
a Floquet-like method (non-Hcrmitian perturbation the- 
ory - NHPt) that treats the xuv field as a monochromatic 
source. This method has been extensively tested in the 
context of x-ray absorption on the Is — > 3» resonance 
of laser-dressed neon, alluded to above @, Ha]- Next, we 
outline a method using direct integration of the time- 
dependent Schrodinger equation (tdse) which achieves 
essentially the same goal using pulses of finite duration. 
The cross section is extracted by projecting out the initial 
state from the final state wave function. The two meth- 
ods, NHPT and the TDSE-projection method, are shown 
to agree when the xuv pulse bandwidth is very small. 
The TDSE-projection approach is, however, potentially 
more flexible in dealing with situations where the dress- 
ing laser couples many states of the atom. We find that 
the TDSE-projection method can be used with a reason- 
able amount of effort to study laser-dressed absorption 
over a wide range of XUV frequencies and dressing wave- 
lengths. As an example of the usefulness of the method, 
we show representative results for several dressing wave- 
lengths between 0.4 and 2.0 /irn. 

We next extend the treatment of the xuv interaction 
to deal with cases where the atomic response varies as a 
function of time or frequency in a non-trivial way. This 
could for instance be because the IR dressing pulse is so 
short that non-adiabatic effects cause the cross section 
to vary substantially over the bandwidth of the pulse. 
Another interesting situation is when the IR pulse is so 
strong that multiphoton processes cause the atom to ex- 
change energy with the light field over a large range of 
frequencies in many different orders of nonlinearity, so 
that it is no longer practical to (artificially) separate the 
linear/non-linear absorption from the driven linear/non- 
linear emission. We therefore develop a time-frequency 
approach to the atomic response, based on the time- 
dependent energy exchange between the atom and the 
light fields. Our method is similar in spirit to the treat- 
ment provided by for instance Tannor [24| and Pollard 
and Mathies [25| with the important difference that we 
are not separating the atomic response into different lin- 
ear and non-linear orders but keeping everything in one 
frequency dependent response function. We find that 
when we use long (~30 fs) xuv pulses we get good agree- 
ment between linear absorption cross sections calculated 
using the time-frequency and the TDSE-projection cross 
section. Having obtained this good agreement over a 
range of frequencies and dressing laser intensities gives us 
confidence that we can calculate the full time-frequency 
response of the dressed atom. 

Finally we show how this time-frequency approach is 
consistent with our solution of the coupled Maxwell wave 
equation (mwe) and the tdse. This allows for a gener- 
alized, ab-initio description of linear and non-linear ab- 
sorption, emission, and phase matching in a macroscopic 
medium. We apply this formalism to studying the prop- 



agation of an xuv pulse in a macroscopic helium gas 
dressed by a moderately intense 11 fs, 1600 nm laser 
pulse. We find that the xuv pulse, which is resonant with 
the Is — > 2p transition in the undressed atom, undergoes 
rich temporal absorption and reshaping dynamics. 

The paper is structured as follows. In Sec. [H] we 
first discuss the three formalisms for calculating absorp- 
tion cross sections of laser-dressed atoms. In Sec. Mil we 
present our framework for the macroscopic calculations. 
Then we use the methods to study laser-dressed helium; 
computational details are given in Sec. IIVI and results 
are presented in Sec. [V] We end the paper with a brief 
conclusion in Sec. I VII 

II. SINGLE ATOM RESPONSE 

This section contains three derivations of the one- 
photon absorption cross section for xuv light interact- 
ing with an atom in the presence of a long wavelength 
dressing field. All three formalisms are based on the sin- 
gle active electron (sae) approximation, and in all cases 
we use linearly polarized fields where the IR and xuv 
polarization vectors are parallel. We use atomic units 
through-out this section [261 ] . 

A. Non-Hermitian Rayleigh-Schrodinger 
perturbation theory 

Our NHPT treatment of dressed xuv absorption is dis- 
cussed in detail in references [H, [III H?} ■ Here we pro- 
vide a brief account to highlight the essential steps in the 
derivation and to facilitate a discussion of the other two 
formalisms. 

In the NHPT formalism, the one-photon xuv absorption 
cross section follows from 



where Jx is the constant xuv photon flux of a continuous 
wave xuv light source [l9j . and Ti is the transition rate 
from the initial state to Rydberg orbitals or the contin- 
uum. The factor of 2 accounts for the number of electrons 
in the atomic orbital which is used as the initial state 1 1 ) . 

To determine T\ with NHPT the full Hamiltonian of an 
atom in two-color light H = Hq + Hi is decomposed into 
a strongly interacting part Hq = Hat + Hem,l + + 
Hem,x that contains the atomic electronic structure if at 
in Hartree-Fock-Slater approximation [28|,[2j|. The inter- 
action with light is expressed in terms of nonrelativistic 
quantum electrodynamics [IH, [l^, [27j ; the free IR laser 
and XUV fields are i?EM,L and i?EM,x, respectively, and 
the interaction of the atomic electrons with the laser field 
is H\ ^. The weak interaction with the xuv light is rep- 
resented by Hi = Hi tX 0. 

Next we represent Hq in a complex-symmetric direct 
product basis of electronic states — without the initial 
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state 1 1) — and photonic number states. In doing so, we 
assume that the initial state and its energy E\ are not 
noticeably influenced by the laser dressing. The matrix 
representation of Hq , 



rr(m)-,(m) 
n C F 



El 



(2) 



is diagonalized, yielding eigenvectors c^ m ' which rep- 
resent the expansion coefficients of new laser-dressed 
states 

|i?(m)) for eigenvalues E^ [H]. Here, m is the 
magnetic quantum number that is conserved for linearly 
polarized light. 

When the Hamiltonian Hq is represented in the new 
basis of laser-dressed states [Eq. fl2J)], the excitation 
or ionization of a ground-state electron of an atom 
due to xuv photoabsorption is described as a reso- 
nance in the spectrum of the non-Hermitian, complex- 
symmetric representation of the Hamiltonian in the basis 
{ 1 1 ) , | ) | VF, m} [Mil . The complex energy of the 
resonance state that 1 1 ) becomes due to the coupling to 
excited states and the continuum via XUV light is usually 
refered to as the Siegert energy [3(1 HI] and satisfies 



E rC s — Er — iTi/2 



(3) 



The real part of the resonance energy is Fr, and 
stands for the transition rate from the ground state to 
a laser-dressed Rydberg orbital or the laser-dressed con- 
tinuum [Eq. @] via photoabsorption. We determine 
the Siegert energy [Eq. ([3])] of the initial state |I) in 
second-order non-Hermitian perturbation theory. The 
total transition rate out of 1 1 ) is given by 



Tt = 2 Im 



(I\Hi |F( ro )) (F^\H! |I) 
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and the absorption cross section is finally obtained from 
Eqs. dH and (jl) as: 



u(u)x) = 87raojx Im 



E 



ha 



F E F 



(m) 
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Here a denotes the fine-structure constant and T> F ' is a 
complex-scaled transition dipole matrix element between 
the initial state |I) and the Fth laser-dressed atomic 
state with projection quantum number m [l9j . 



B. Projection treatment of xuv absorption 

As an alternative to the treatment above, we can ob- 
tain the linear absorption cross section by a direct so- 
lution of the tdse in the SAE approximation [34j . The 
cross section is extracted by projecting the final state 
wave function obtained at the end of a finite pulse onto 
the initial wave function. As such, we avoid calculating 
the dressed states directly, making explicit use of only 
the laser-free initial and final states. 



To simplify the treatment of finite duration pulses 
when using the projection method, we replace the 
quantum electrodynamic treatment of xuv radiation in 
Sec. Ill Al by a semiclassical description of light |35j . We 
begin by choosing the vector potential of the xuv light 
of carrier frequency u>x to be 



£x(t) 



sin(w x t) e x 



(6) 



The electric field of the xuv light field is then given by 
derivative with respect to time, £x(t) — —dAx{t)/dt: 



£x(t) 



£x(t) cos(u;x t) 



i ggxft) 
dt 



sin(wx t) 



(7) 

Here, £x(t) = y/8nalx(i) is the envelope of the xuv 
pulse and Ix(t) is its cycle- averaged intensity. Our spec- 
ification of the vector potential in Eq. [5] ensures that the 
integrated electric field and the vector potential at the 
end of the pulse, A(tf), are zero when £x(t) is zero at 
the initial and final times. It leads to the second term on 
the right hand side of Eq. [7] which is a small correction 
of order Awx/wx near the center of the pulse, for pulses 
with a bandwidth of Aujx- By ensuring that A(tf) = 
we obtain results which are independent of the electro- 
magnetic gauge. 

For a Gaussian envelope pulse with a full width at 
half maximum (fwhm) duration of tx, the bandwidth of 
the pulse is given by Aujx — 41n2/rx- In our calcula- 
tions, we first specify Acjx and this dictates the value of 
tx- The Gaussian envelope is then approximated by a 
trigonometric pulse (3r| : 

Ix(t) = Ix,o cos 2n (g) - |i|) = J X ,0 9n(t) , (8) 

with an integer n > and the Heaviside 9 function [37| . 
The total pulse duration is defined as 



T„ 



TTTx 



2 arccos 2 ^ 



(9) 



The envelope ((8]) converges rapidly to a Gaussian func- 
tion in the limit lim g n (t) — > exp ( — 41n2( — ) 2 ) . Using 

the approximative function JSJ instead of a true Gaus- 
sian function has the advantage that it goes to zero on a 
finite support, which allows us to satisfy the requirement 
A(tf) = exactly. 

In the TDSE-projection formalism we also need to dress 
the atom with a laser field with frequency ul . We do this 
by using a laser field of the form 



£ L (i) = £h(t) e L sin(wLt) 



(10) 



The envelope function £h(t) is now a trapezoidal pulse 
with a linear ramp of one optical cycle at each end and 
a flat section that completely spans over the XUV pulse. 
The pulse contains an integer number of laser cycles, so 
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we again obtain zero vector potential at the end of the 
dressing pulse. We assume this field is too weak to excite 
or ionize the atom on its own, an assumption that we can 
explicitly check by running the calculation once without 
the xuv field. 

To calculate the cross section for absorption we begin 
with the atom in its ground state at time to and 
use the grid-based methods of reference [HI to propagate 
the wave function forward in time until the end of the 
combined xuv and dressing pulse at time tf. At this time 
we calculate the probability that the atom has remained 
in its ground state Pi(tf) by projecting the final wave 
packet \ip{tf)) onto the initial wave packet: 



(ii) 



Given Pi(oo) = P\{tf) from the tdse calculation, we 
obtain the probability that the atom is excited or ion- 
ized from 1 — Pi(tf). Because we are dealing with a 
one-photon absorption process where we assume the in- 
tensity is well below saturation, a linear relation be- 
tween xuv absorption rate and xuv photon flux holds: 
r(i) = a(u)x) Jx{wx,t). We use this assumption to 
transform the probability to absorb an xuv photon into 
an expression for the cross section: 



2[1 



Pi(oo)] = a{u> x ) J Jx{u>x,t) dt 



(12) 



This is equivalent to the steady state expression in 
Eq. (fTJ : The factor 2 again stems from the two electrons 
in the spatial orbital I which contribute equally. 

The underlying assumption in Eq. (I12[) is that we can 
calculate the absorption cross section <r(u>x) for a small 
range of frequencies Awx around wx by calculating the 
response of the atom to a pulse of bandwidth Awx . For 
Eq. (fT2"j) to be meaningful, the cross section needs to be 
approximately constant over the bandwidth wx of the 
pulse. For a low-bandwidth pulse, we can further use the 
relation Ix(t) ~ wx Jx(wx,t) between photon flux and 
intensity. Then, the time-integral on the left-hand side 
can be solved analytically for the pulse shape ([7]). In this 
way, we find the xuv absorption cross sections <j(ojx) 
from Eq. (|12|) by dividing the probability to excite an 
atom out of the ground state 1 — P\{tj) by the integral 
over the xuv flux. 

As we stated in the introduction, though we expect 
that the two methods for calculating the frequency- 
dependent absorption cross section should agree, the 
tdse projection approach is potentially more flexible in 
dealing with situations where the dressing laser couples 
many states of the atom which forces the Hamiltonian 
matrix in Eq. [5] to be very large. 



C. Time- frequency treatment of ultrafast xuv 
absorption 

In this section we extend the treatment of the xuv in- 
teraction to deal with cases where the atomic response 
varies as a function of time or frequency in a non-trivial 
way. This could be because the dressing IR pulse is 
substantially shorter than the xuv pulse, or when non- 
linear interactions would cause the atom to exchange en- 
ergy with multiple xuv frequency components in differ- 
ent non-linear orders. 

We start by deriving a frequency-dependent response 
function S(u>) from the time-dependent energy exchange 
between the atom and the light field. S(lj) is defined so 
that when integrated over all frequencies, it yields the 
total excitation probability. This includes excitation to 
continuum states, i.e., ionization. We can then express 
the total energy gained by the atom from the light fields, 
AE, as the sum over the frequency dependent excitation 
probability S(ui) times the photon energy: 



AE = 



lo S(lo) dio. 



(13) 



To calculate the response function we use that the total 
atomic energy gain can also be expressed as a sum over 
the rate at which energy is gained: 



AE = 



lo S(lo) dio 



dE 
~dt 



dt. 



(14) 



We calculate this rate directly from our one electron 
Hamiltonian, H — Ha + £{t) z, as: 

f-lww-wigw-wf. («) 

We note that £ (t) is the full electric field consisting of 
the sum of the dressing laser and the xuv fields. This 
means that we are simultaneously treating the exchange 
of energy between the atom and all frequencies of the 
light field. In the following we will denote (z)(t) by z(t). 
The time-dependent dipole moment is related to z(t) by 
d(t) = —z(t) for a single electron. We now calculate AE: 



AE = 



z(t)-dt 



(16) 



= -J w2Im{,?(u>)£*(u;)} dw . (17) 

In this derivation we have used that both z{t) and £{t) 
are real functions of time so that z(— lo) = z*(uj) and 
£{— lo) — £*{w). Using Eq. (fT4"]) we then have an expres- 
sion for the response function: 

S+(lo) = -2Im{z(w)£*(w)j lo > 0, (18) 

where the -I- subscript on S + (lo) explicitly indicates that 
we are only integrating over positive frequencies. 



5 



We calculate the dipole spectrum in the SAE approx- 
imation dsAE{^>) via the time-dependent acceleration 
a(t): 



d 2 z 



-mw,[H,z\}\m), (19) 



The dipole spectrum is then given by dsAE{w) — 
d(uj)/uj 2 , where a(w) denotes the Fourier transform of 
a(t). The full (two-electron) dipole moment is d(uj) — 
2d sae (w). 

In the weak-lR limit where it is meaningful to talk 
about an absorption cross section, we can write the 
frequency-dependent energy exchange function ujS(oj) by 
means of a generalized cross section a(ui) and the spectral 
energy density of the electric field, luJ(uj). The spectral 
flux J(w) is defined as [53j : 



J(w) = 



1 



Anauj 



£(uj) 



(20) 



This means that once we calculate the response function 
S(lu), the generalized cross section is given by: 



<r(w) 



(21) 



Inserting the response function from Eq. (|18|) we obtain 
the cross section, now defined for both positive and neg- 
ative frequency components: 



<t(uj) = 8nauj Im 



\ d~SAE{u) 

\ £{u) 



(22) 



This equation is the generalized, time-frequency, multi- 
mode equivalent of Eq. ((TJ) which was derived for the 
steady-state case. 

To calculate the generalized cross section in Eq. (I2"2j) , 
and the macroscopic polarization field described in the 
following sub-section, we multiply the time-dependent 
acceleration in Eq. (|19l) with a window function W(t), 
aw(t) — a(t)W(t) and calculate dsAE^) from the 
Fourier transform of aw if)- In Eq. (|2"2"j) we also cal- 
culate £{ui) from W{t)£{t) for normalization purposes. 
The window function on the time-dependent accelera- 
tion is necessary in particular in those cases where the 
xuv light is resonant with an atomic transition. The 
xuv light then induces a strong coherence between the 
ground state and the excited state which in the numer- 
ical calculation will go on "ringing" until long after the 
xuv pulse is over. This ringing does not correspond to 
stimulated emission or absorption of xuv radiation. The 
window function we use is a trigonometric function as 
given in Eq. © and is in general chosen to have the same 
fwhm duration as the longer of the IR and xuv pulses. 
The choice of window function has some influence on the 
value of the cross section for the un-dressed atom around 
the field-free resonances. When the atom is laser-dressed 



so that the xuv light is no longer absorbed as strongly, 
the ringing is strongly suppressed by the laser field and 
the influence of the window function is very small. 

It is interesting to note here that for intense or few- 
cycle IR fields, and/or for multi-mode xuv fields, the sign 
of the response function S+(w) (and therefore the sign of 
the generalized cross section) for a particular frequency 
u> in Eq. (fT51) can be positive or negative. When S+(w) 
is positive the atom will predominantly absorb light of 
that frequency, and when S+(lj) is negative the atom will 
predominantly emit light of that frequency. This makes 
the response function a powerful tool for studying the 
dynamics of the light-atom energy exchange, in particular 
in combination with a sliding time-window on the time- 
dependent acceleration. This would in principle allow 
for the time-resolution of when different frequencies are 
absorbed or emitted during a dynamical process. We will 
discuss a simple application of this in connection with the 
macroscopic reshaping of an xuv pulse presented in the 
Results section. 



III. MACROSCOPIC RESPONSE, INCLUDING 
ABSORPTION 

As we will show at the end of this section, the relation- 
ship derived in the previous section, between the dipole 
spectrum driven by an arbitrary pulse and the absorp- 
tion cross section for the frequencies contained in that 
pulse, is consistent with our general framework for the 
interaction between an ultrafast, multi-color pulse and 
a macroscopic medium. This framework consists of the 
coupled solutions of the MWE and the tdse for all fre- 
quencies u> of the electric field E(lu) of the multi-color 
pulse. We will express all quantities in SI units in this 
section. In a frame that moves at the speed of light, and 
in the slowly evolving wave approximation which works 
well even for few- femtosecond pulses [Hj , the mwe takes 
the following form: 



V 2 ± £(u>) 



2icu d£(ui) w 



dz 



-rfPH+P^M). (23) 



The electric field £{uS) and the source terms P(co) and 
Pion{u) are also functions of the cylindrical coordinates r 
and z. We solve this equation by space- marching through 
the helium gas, at each plane z in the propagation direc- 
tion calculating the response terms P(ui) and Pi on (u;) via 
numerical integration of the tdse, and then using them 
to propagate to the next plane in z. The macroscopic 
polarization field P(u) is calculated from two times the 
one-electron single atom dipole moment dsAE^)- 



P(w) = 2pd SAE {w) 



2pe 



a{t)W{t)e %ut At, 



where p is the atomic density, a(t) is the time-dependent 
acceleration calculated as described in Sec. Ill CI As the 
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driving field for the tdse calculation we use the evolving 
electric field £ (t) at the plane z. This means that P{oj) in 
general includes both the linear and nonlinear response 
of the atom to the multicolor field. The term Pi on {ui) is 
due to the space- and time-dependent free-electron con- 
tribution to the refractive index and is also calculated 
within the SAE-tdse, see [39|. This term is very small in 
the cases considered in this paper and we will ignore it 
hereafter. 

By calculating the source terms in each z-plane and 
using them to propagate to the next z-plane, we are cou- 
pling both the linear and non-linear response generated 
in one step back into the full electric field so that it can 
contribute to the driving electric field in the next step. In 
much of the work described in the literature, see for in- 
stance [38U42| , the non-linear response is separated from 
the linear response, and the propagation of the newly 
generated radiation (via nonlinear processes) is separated 
from the propagation of the driving field. Absorption 
and dispersion of different frequency components of the 
light fields are then added separately, typically using tab- 
ulated, frequency-dependent values. It has been shown 
in a number of papers that such an approach offers a 
very complete description of both the generation of new 
frequencies via nonlinear processes, and the macroscopic 
effects of phase matching and ionization-driven reshap- 
ing of the ultrafast propagating pulse [38l - l42j . However, 
it cannot describe ultrafast or dynamical reshaping of 
the xuv pulses driven by for instance absorption, dis- 
persion, or laser-induced transparency. More generally, 
processes that are due to the combined response to the 
strong dressing or driving laser field and the weaker xuv 
fields are not described in a self-consistent manner be- 
cause the generated radiation is not included into the 
driving field. 

In the following we will argue that the approach pre- 
sented in this paper, which allows us to calculate the 
non-linear response of the dressed atom, also allows us 
to describe the absorption and dispersion of the ultra- 
fast pulses in a self-consistent manner, to within the SAE 
approximation. Let us first rewrite the macroscopic po- 
larization field as: 

P{u) = pd{uj) = P [Re(§^-) + flm(f|)]^H (25) 

The last term on the right is proportional to the general- 
ized cross section in Eq. (|22[) . By inserting this expression 
into the MWE in Eq. (J22J) we get: 



d£{co) 
dz 



-^pRe(M)f( w )-ipcr( W )f( W ). 
£oc £{oj) 

The second term on the right hand side clearly will lead 
to absorption at frequency ui with absorption coefficient 
pd{ui) when <t(cj) is positive, which is the case in the 
weak field limit when the atomic response is linear. In 
this linear case, the first term on the right hand side 
can likewise be interpreted as a generalized expression 



for the dispersion experienced in the gas medium, with 
the frequency dependent correction to the refractive in- 
dex given by Ari(w) = ^R-Kf^ )■ The strength of our 
time-dependent approach is that even when the driving 
field is strong enough to induce non-linear processes, we 
are able to treat all of the linear and non-linear processes 
within one time-dependent calculation, rather than ar- 
tificially separating processes of different non-linearities 
and assigning them a frequency- and intensity-dependent 
weight. 

IV. COMPUTATIONAL DETAILS 

Computations with the time-independent theory of 
Sec. Ill Al were carried out with the DREYD computer pro- 
gram from the fella suite [43j. The computational pa- 
rameters are specified in analogy to Rcf. [19(. However, 
in this work, we do not rely on the Hartree-Fock-Slater 
mean-field approximation [28, 29] to describe the atomic 
electronic structure. Instead, we use a pseudopotential 
for helium, constructed from the ground state Hartree- 
Fock potential, calculated on a very fine radial grid by 
standard iterative methods 44j . We set the K edge of 
helium to the value of E ls = -24.5786 eV. Next, the 
radial Schrodinger equation is solved with the pseudopo- 
tential where the solution, the radial part of the atomic 
orbitals, is represented on a grid with a radius of 60 ao us- 
ing 3001 finite-clement functions. From its eigenfunctions 
we choose, for each orbital angular momentum I, the 
100 functions which are lowest in energy to form atomic 
orbitals [19j. In doing so, we consider spherical harmonics 
with up to I = 7 [37], Ell . Continuum electrons are treated 
with a smooth exterior complex scaling complex absorb- 
ing potential 46-48] which is parametrized with the com- 
plex scaling angle 9 = 0.13 rad, a smoothness of the path 
of A = 5 <2q 1 , and an exteriority of ro = 10 ao [l9j ]. There 
is only radiative decay of singly excited states of helium 
with comparatively long lifetimes to all other time scales 
in the problem; therefore, we set the linewidth of a if va- 
cancy in helium to zero. Finally, we diagonalize the in- 
volved Floquet-type matrices to obtain the cross section. 
Without the laser field this is done exactly; when the 
laser is present we use 4000 Lanczos iterations [l8j |. 

Computations with the TDSE-projection method of 
Sec. Ill Bl were carried out with a one-electron tdse 
solver code which is based on the algorithms described 
in Ref. (34)]. The same potential used above is trans- 
ferred to a radial grid with spacing of 0.2 ao and used for 
the TDSE-projection and fully time-dependent computa- 
tions (see below). The interpolation of the Hartree-Fock 
potential onto the coarse grid used for the TDSE propa- 
gation introduces a small error in the helium Is ionization 
potential which we correct by slightly changing the po- 
tential at the first grid point [3J]. The pulse shape is 
given by Eq. ([7][9|) with n = 6. This means that the to- 
tal propagation time is 4.67tx- Typically we use a box 
of 200 au in size, with a 50 au absorbing boundary at 



7 



the outer edge 34]. The maximum angular momentum 
and time step size are adjusted to achieve convergence. 
We use £ max = 8 and 1500 steps per dressing laser cycle. 
In some cases where the xuv pulse was very long or the 
xuv wavelength was very close to the ionization thresh- 
old, the box size was increased to 1000 au ensure that no 
wave function amplitude that might reflect from the ab- 
sorbing boundary could interfere with amplitude excited 
at a later time. We specify the bandwidth Awx of the 
XUV pulse instead of the FWHM duration tx as in Eq. (jHJ . 
For our (approximately) Gaussian pulse, we use the time- 
bandwidth product rx Acjx = 4 In 2 to convert between 
the two quantities [49| with Awx = 0.05 eV which corre- 
sponds to a duration of rx = 36.5 fs. We investigated the 
dependence of the absorption cross sections on the inten- 
sity of the XUV light; to a very good approximation, we 
find a linear relationship as should hold for a one-photon 
absorption process ([T]). An xuv intensity of 10 10 is 
employed in Figs. 1-4. 

The calculations with the time-frequency method of 
Sec. Ill Cl were performed with the TDSE solver described 
above. For the calculations in Figs. 1(b) and (c) and the 
inset in Fig. 2 we have used t max = 8 and approximately 
4000 steps per cycle of the dressing laser field. The size 
of the radial grid was 150ao (using 750 points) with a 250 
point absorbing boundary. The intensity envelope of the 
m pulse is cos 4 (/3£/i7_r), where tir is the fwhm duration 
of the IR pulse and /3 = 2 arccos(0.5 1 ^ 4 ). The intensity 
envelope of the xuv is usually chosen to be the fourth 
power of the IR envelope (to be consistent with the xuv 
being a high order harmonic produced by the IR pulse). 
This gives a fwhm pulse duration for the xuv pulse of 
about half that of the IR pulse. The window function 
discussed in Sec. Ill CI is a Hann window with a fwhm 
duration very close to that of the IR pulse. The window 
function was chosen such that the long-pulse calculation 
in Fig. 1(b) can be compared to those in Fig. 1(a): the 
fwhm bandwidth of the windowed acceleration spectrum 
aw{u) has the same 0.05 eV bandwidth as the TDSE- 
projection approach. 

For the mwe-tdse calculations in Fig. 6 we employ two 
time scales. One time scale defines the spectral resolu- 
tion of the macroscopic, propagating, electric fields. This 
time scale typically extends to ±3 times the fwhm of the 
longest of the IR and xuv pulses and contains approxi- 
mately 5500 time points. The other time scale is used 
for the TDSE solution and extends only over the finite 
duration of the longest pulse, and typically contains 6000 
points per IR laser cycle. The macroscopic length scales 
cover 160 /im in the radial direction, with 200 grid points, 
40 of which contain an absorber that prevents reflections 
from the edge of the grid, and 1 mm in the propagation 
direction, with 600 grid points. In the propagation direc- 
tion we only evaluate the dipole moment every 20 steps, 
and rescale the response to the appropriate density and 
phase in between, see [3!| for details. The initial spatial 
distribution of both the xuv and the dressing laser beam 
is Gaussian. The xuv beam has a confocal parameter of 



10 cm and a corresponding focal diameter of 60 /im. The 
1600 nm dressing pulse has a confocal parameter of 2 cm 
and a focal diameter of 140 /im. This means that in the 
spatial dimension, the xuv beam is always overlapped 
with the IR beam. The IR beam changes only marginally 
during the propagation in the helium gas. 



V. RESULTS AND DISCUSSION 
A. Single atom absorption cross sections 

The helium absorption cross section for linearly po- 
larized XUV light, in the absence of laser light, is dis- 
played in Fig. [TJ In part (a), we compare results from 
dreyd [43| with results of the TDSE-projection method 
[Eq. (|12l) ]. To be able to compare these two results, we 
have convoluted the DREYD cross sections with a Gaus- 
sian with the same bandwidth of Awx = 0.05 eV that 
was used in the TDSE calculation. This leads to good 
agreement between the two results. We note that the 
presence of a spectral bandwidth in both calculations 
means that we are only able to resolve spectral features 
to within 0.05 eV. The peaks at 21.1068 eV, 23.0416 eV, 
23.7162 eV, 24.0273 eV, . . . stem from Is 2 -> lsnp tran- 
sitions with n £ {2,3,4,5, . . .}. In parts (b) and (c) we 
show cross sections calculated using the time-frequency 
approach leading to Eq. (|2"2"j) . around the 2p and the 6p 
and 7p states. These calculations were done using an ex- 
tremely weak 764 nm IR pulse and harmonics 13 (b) or 
15 (c) of the IR frequency 54j. Harmonic 13 is resonant 
with the 2p state and harmonic 15 is in between the 6p 
and the 7p states. We show the results of using three dif- 
ferent xuv pulse durations (30 fs, 15 fs, and 7.5 fs). The 
IR pulse has twice the duration of the xuv pulse and an 
intensity of 10 8 (low enough that it does not influ- 
ence the cross sections). The 30 fs calculation leads to a 
0.05 eV bandwidth of the dipole moment around the 2p 
state, after applying the time-domain window function 
discussed in Sec. [TV] The calculated cross section is in 
reasonably good agreement with the results in (a). The 
shorter xuv pulses lead to broader absorption cross sec- 
tions. For the 15 fs xuv pulse the 6p and 7p states can 
still be distinguished as separate features in the absorp- 
tion spectrum. Using 7.5 fs XUV pulse the cross section 
can be calculated over a much larger frequency range, 
spanning both below and above the ionization threshold, 
and as a consequence one can no longer distinguish the 
6p and 7p states. The value of the cross section in this 
calculation is in good agreement with the value in Hcnke 
et al. [50( of 7.5 Mbarn just above threshold, as we expect 
when using pulses that span the ionization threshold. 

In Fig. [2] we show how the xuv cross section changes 
when the helium atom is exposed to an infrared laser field 
with an intensity of 10 12 and a wavelength of ~ 800 
nm. The main figure again compares the results from 
dreyd [43] and the TDSE-projection method (which have 
again both been calculated/convoluted with a 0.05 eV 
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FIG. 1: (Color online) The xuv absorption cross section of 
a helium atom, (a) The dashed red lines were obtained with 
dreyd [43] the solid black lines were obtained using the TDSE- 
projection method, (b) and (c) show close-ups of the cross 
section around the 2p (b), 6p, and 7p (c) states calculated 
with the time-frequency method [Eq. (I22| for different xuv 
pulse durations. The results obtained using 30 fs, 15 fs, and 
7.5 fs xuvpulses are shown in black (circles), red (squares), 
and green (open diamonds), respectively. 



I 6- 
S - 

o 5- 
o 

a 4- 
c 

%*- 

CD 

s 2- 

P 1 - 



o 



_ 

20 



21 



A i 




22 23 
Photon energy (eV) 



Photon energy (eV 



24 



25 



FIG. 2: (Color online) The xuv absorption cross section of 
a helium atom dressed by an intense 800 nm ir laser pulse 
with a peak intensity of 10 12 The dashed red lines were 
obtained with dreyd [43J; the cross section obtained from 
the TDSE-projection is plotted with dashed black lines. The 
inset shows the cross section calculated using Eq. (|22[> . using 
an xuv pulse duration of 7.5 fs (xuv pulse duration 7.5 fs), 
and an 764 nm ir pulse with a duration of 15 fs and a peak 
intensity of 10 s ^ (solid black curve), 10 11 ^ (dashed red), 
5 x 10 11 (dot-dashed green), or 10 12 ^ (dotted blue), 
respectively. 
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bandwidth) , and they are found to be in good agreement 
over a broad energy range. 

The absorption of the dressed atom in Fig. [2] changes 
significantly from the undressed case, although many of 
the field-free resonances can still be recognized. The 
peak due to the 2p state has broadened and shifted to 
lower energy, whereas the higher np peaks are shifted to 
higher energies. In addition, several new absorbing fea- 
tures have appeared between 21 eV and 22 eV. The inset 
in Fig. [5] shows cross sections for the dressed helium atom 
calculated using the time- frequency approach of Eq. (f22j) . 
for an xuv pulse duration of 7.5 fs. The 764 nm ir pulse 
duration is 15 fs and the ir peak intensity varies between 
10 s (undressed, as shown in Fig.rjjb)) and 10 12 
The inset details the shift and broadening of the 2p res- 
onance as the dressing laser intensity is increased. We 
have chosen to use the 7.5 fs xuv pulses for these cal- 
culations in order to be able to cover the shift of the 2p 
resonance within the bandwidth that can be addressed 
within Eq. ff22]>. 

The extra peaks in the 800 nm dressed-atom cross sec- 
tion shown in Fig. [5] result from complex multiphoton 
effects and do not have a straight forward interpreta- 
tion. Other dressing-laser wavelengths offer more insight 
into the non-linear optics driven by the two-color field. 
We first show two figures exploring the impact of the 
dressing laser wavelength on the xuv absorption cross 
section. The wavelengths we have used are listed in Ta- 
ble U together with the corresponding photon energies. 
All of these wavelengths can be produced from standard 



TABLE I: Correspondance between wavelenght and photon 
energy for the involved laser light. 



Ti:Sapphire high power, short pulse laser systems via fre- 
quency mixing in nonlinear materials. 

The laser-dressed xuv absorption cross section, calcu- 
lated using the TDSE-projection method are displayed 
in Figs. [2] [3j and [4] for several laser wavelengths at an 
intensity of 10 12 For 400 nm and 500 nm light, we 
see only a moderate impact of the laser dressing. The 
impact is mostly on the 2p state, as the largest dipole 
coupling exists to other close-by Rydberg states. 620 nm 
and 800 nm light exhibit complex multiphoton effects 
which manifest in complicated multi-peak structures in 
the cross sections. The dressing pulses with longer wave- 
lengths all induce systematic behavior. In all three cases, 
the single Is 2 — > ls2p transition in Fig. Ufa) (without 
dressing) is split into two lines by the laser in Fig. [3J 
the transitions from the Is orbital into higher Rydberg 
orbitals are replaced by a continuous, weak absorption 
feature. 

We would like to elucidate the origin of the double- 
peak feature around 21 eV in the long- wavelength series 
shown in Fig. [4] It is much simpler than the correspond- 
ing feature for the wavelengths in Figs. [5] and [3] To this 
end, we make a A-type model for helium which is shown 
in Fig. [5] It comprises the ground state of helium and 
the ls^ 1 2p, Is 2s excited states. The laser photon 
energy is denoted by ojl whereas r ls -i 2s and r ls -i 2 p 
are parameters for the laser-induced decay widths of the 
respective excited states. The overall agreement of the 
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FIG. 3: (Color online) Laser-dressed xuv absorption cross 
section of helium for 400, 500, and 620 nm laser wavelengths 
at a laser intensity of 10 12 -^r. 
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FIG. 4: (Color online) Laser-dressed xuv absorption cross 
section of helium for 1400, 1600, and 2000 nm laser wave- 
lengths at a laser intensity of 10 12 (red solid curves). 
The results of the simple three-level model of Fig. [S] are in- 
dicated by the dashed black curves for r ls -i 2p = 0.1 eV and 
r ls -i 2s =0.05eV. 



model curves with the ab initio data in Fig. 2] is sat- 
isfactory. The reason for the success of the three-level 
model is — as in Ref. 0] — the fact that the splitting be- 
tween the 2s and 2p Rydberg orbitals in helium is 0.84 eV, 
i.e., the laser is almost in resonance with this transi- 
tion, within the laser-induced line widths, for midinfrared 
wavelengths [Table []. Furthermore, the other levels of 
helium couple only weakly. 

The A-type model explains the double peaked struc- 
ture in Fig. 0]in terms of a splitting of the Is -1 2p and 
Is -1 2s states into an Autler-Townes doublet. This fea- 



He 2s 



FIG. 5: (Color online) A-type three-level model for helium. 
The laser photon energy is cjl and the xuv photon energy 
is cjx- The laser-induced decay widths of the Is -1 2p and 
Is -1 2s excited states are denoted by r ls -i 2p and r ls -i 2s , 
respectively. 



ture raises the possibility that the dressing laser could be 
used to induce transparency to the xuv radiation tuned 
to the Is 2 — > ls2p transition. A similar mechanism was 
found for the suppression of resonant absorption of x rays 
in neon ll5T - fl7] ] ■ argon [TH, and krypton [19[ atoms and 
called eit for x rays Q . In the next section we study the 
analogous effect in helium. 



B. xuv pulse shaping in a macroscopic medium 

In this section we present an application of the time- 
frequency approach to absorption in a macroscopic non- 
linear medium. We study how the laser induced trans- 
parency discussed above may be used to temporally con- 
trol the xuv pulse shape in a helium gas, in analogy with 
the EIT for x-rays discussed in In that x-ray study, 
the absorption was exclusively described in terms of an 
intensity-dependent absorption cross section, which then 
in turn enforces a one-to-one mapping of the absorption 
to time through the intensity. This does not allow for 
truly dynamical effects. In addition, the intensities we 
explore here are much lower that those used in the x-ray 
study, which means that ionization of the Rydberg states 
does not play a large role. 

We calculate the electric field of a combined two-color 
xuv-ir pulse after propagation through a 1 mm long he- 
lium gas jet with a density of 1.5 x 10 17 cm -3 (6 mbar 
at room temperature). We solve the coupled MWE-TDSE 
in the form of Eq. (j2U)l as described in section IIII1 see 
also [1^]. The initial xuv pulse has a wavelength of 
58.7 nm (21.1 eV, resonant with the Is 2 — > ls2p tran- 
sition), a pulse duration of 61 fs and a peak intensity of 
10 7 W/cm 2 . The 1600 nm dressing pulse has a peak in- 
tensity of 10 12 W/cm 2 . We have used different IR pulse 
durations between 122 fs and 11 fs. We have checked 
that the reshaping discussed below is no different when 
we use a higher xuv intensity of 10 10 W/cm 2 . 

Fig. [5] shows the radially integrated spectrum (a) and 
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FIG. 6: (Color online), m-assisted xuv absorption in 1 mm 
long macroscopic helium gas with a density of 4 x 10 16 cm -3 . 
The initial 61 fs xuv pulse is resonant with the 2p state of 
the undressed helium atom. We show the xuv spectrum in 
(a) and time profile in (b), both before (solid lines) and af- 
ter propagation. Final profiles at the end of the un-dressed 
medium are shown with dotted red lines, final profiles at the 
end of the medium dressed by an 11 fs, 1600 nm ir pulse with 
a peak intensity of 10 12 W/cm 2 are shown with dashed blue 
lines. In (b) we also show the final time profile when the in- 
tensity of the dressing pulse is 2 x 10 12 W/cm 2 (thin green 
line) . 



time profile (b) of the xuv pulse before and after prop- 
agation through the helium gas. When the atoms are 
undressed the xuv radiation is strongly depleted via the 
resonant absorption, as is shown by the dotted red lines 
f55j . The absorption length at this atomic density is less 
than 0.1 mm. During the first few absorption lengths 
the xuv yield decreases exponentially. The large disper- 
sion across the resonance, and to a lesser extent the fre- 
quency dependence of the absorption cross section, sub- 
sequently leads to reshaping of the depleted beam upon 
further propagation in the gas. This causes the double- 
peaked shape of the spectrum emerging at the end of the 
medium. The time profile of the final XUV field is cor- 
respondingly irregular, as seen by the red dotted line in 
Fig.lHfb). 



We then apply a 1600 nm, 10 12 W/cm 2 dressing pulse 
which is much longer than the xuv pulse (123 fs vs 61 fs). 
This means that the xuv pulse encounters a sample of 
strongly dressed atoms which are no longer resonant with 
the xuv energy, see Fig. 4, and the gas is therefore trans- 
parent to the xuv light. The spectrum of the final XUV 
pulse is nearly indistinguishable from the initial spectrum 
and is not shown in Fig.^a). The final xuv pulse shape 
is also nearly identical to the initial pulse shape except 
for a 1.6 fs delay caused by the different group veloci- 
ties of the IR and the xuv pulses (also not shown in the 
figure) . 

Next, we apply an 11 fs dressing IR pulse which is sub- 
stantially shorter than the xuv pulse. This means that 
the dressing pulse turns on and off within the FWHM du- 
ration of the xuv pulse, thereby strongly coupling the 
2s and 2p states in a dynamical manner. The final spec- 
tral and temporal xuv profiles at the end of the medium 
are shown with dashed blue lines in Fig. EJa) and (b). 
The xuv time profile at the end of the medium is dom- 
inated by an approximately 10 fs pulse, superimposed 
on a much weaker longer pulse, and the corresponding 
xuv spectrum has broad shoulders at frequencies sub- 
stantially beyond the initial xuv bandwidth. We note 
that the final xuv pulse is not symmetric around time 
zero and that in particular, the short sub-pulse is delayed 
by approximately 4 fs from the center of the dressing 
IR pulse. We attribute this to the complicated absorp- 
tion and emission dynamics driven by the two-color pulse 
as explored in more detail next. 

The time-dependent acceleration driven by the initial 
two-color pulse is shown in Fig. solid black line. We 
are showing the envelope of the acceleration to avoid the 
fast oscillations at the resonance frequency. On the ris- 
ing edge of the xuv pulse, before the ir pulse turns on, 
this acceleration represents the absorption of the xuv 
light via population transfer to the 2p state. This means 
that when the IR pulse arrives there is already popula- 
tion in the 2p state which will couple strongly to the 2s 
state. This causes the suppression of the acceleration 
which starts around t = —5 fs. The acceleration then 
has a revival around 4 fs before it is suppressed again. 
By calculating the generalized cross sections of Eq. (|22l) 
separately for the peaks around t = —9 fs and around 
t = +4 fs, we find that they have opposite signs. This 
means that whereas the dipole response early in the pulse 
and up until t m — 5 fs is due to the absorption of the xuv 
light, the peak in Fig. [7] around t — +4 fs corresponds 
to emission of xuv radiation. We interpret this behavior 
as coming from Rabi-like oscillations of the excited state 
population between the 2s and the 2p states, driven by 
the IR field. The emission happens when the population 
returns to the 2p state while the strong ir field is still 
on. This interpretation would predict that at higher IR 
intensity the Rabi cycling should be faster. We indeed 
find that if we increase the IR intensity to 3x 10 12 W/cm 2 , 
the time-dependent acceleration has two revivals within 
the IR pulse duration, both corresponding to emission 
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FIG. 7: (Color online). Single atom time-dependent accelera- 
tion driven by the initial xuv-m pulse in Fig. [5] The duration 
of the two pulses are 61 fs and 11 fs, respectively, and the xuv 
intensity is 10 7 W/cm 2 . In solid black line we show the result 
of using an ir intensity of 1 x 10 12 W/cm 2 , in dashed red line 
the ir intensity is 3 x 10 12 W/cm 2 . 

(dashed red line in Fig. [7]). The Rabi oscillation period 
for resonant population transfer between the 2s and the 
2p states is approximately 10 fs (6 fs) for a constant in- 
tensity of 1 x 10 12 W/cm 2 (3 x 10 12 W/cm 2 ). This is in 
good agreement with the time scale of the oscillations in 
the acceleration seen in Fig. 7, especially considering that 
the IR intensity is changing rapidly between t = — 10 fs 
and t = 10 fs. 

Finally, returning to Fig. ^[b) and the xuv pulse that 
emerges from the helium gas dressed by the 11 fs IR pulse, 
we can now attribute the delay of the short xuv pulse to 
the excited state dynamics in the strongly dressed atomic 
gas. Fig. [5] shows the evolution of the xuv time profile 
shown in Fig. [51 as a function of propagation distance. 
The complicated atomic response shown in Fig. which 
includes absorption at early times and emission around 
t = 4 fs, is reflected in the propagating XUV electric field. 
After the first few hundred microns of propagation the 
xuv time profile has been substantially depleted on the 
rising edge and is dominated by a much shorter pulse 
peaking shortly after t = 0. We note that this in turn will 
change the atomic response from that plotted in Fig. [7] 
since the dressing IR pulse and the xuv pulse are then 
more comparable in duration. 



VI. CONCLUSION 

In this paper, we have investigated the response of 
laser-dressed helium atoms to xuv radiation, within the 
SAE approximation. In particular, we have focused on 
the calculation of absorption cross sections and their ap- 
plication to absorption in a macroscopic medium. 

First, we introduced a time-independent method based 
on NHPT. The interaction with the XUV light was treated 
in terms of a one-photon process while a Floquet-like ap- 
proximation was used to describe the impact of the dress- 
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FIG. 8: (Color online). Evolution of xuv time profile during 
propagation through the macroscopic helium gas. The time 
profiles at different propagation distances have been displaced 
vertically, starting with z = at the bottom to z = 1.0 mm 
at the top, in increments of Az = 0.01 mm. 



ing laser. Second, we devised a time-dependent method 
to compute the cross section with using a direct inte- 
gration of the tdse and projection of the final wave 
packet onto the initial atomic state. We showed that 
the projection-based approach, which was implemented 
using finite pulses, yields the same results compared with 
the time-independent results. Third, we presented a ver- 
satile time-frequency approach to evaluating an atomic 
response function which can be used even when the dress- 
ing laser pulse is so short that it introduces transient ef- 
fects, or in cases where the atom exchanges energy with 
multiple frequency components of the multi-color light 
field. We showed that this method, when used to calcu- 
late linear absorption cross sections, agrees with the first 
two. Finally, we showed that this third approach can be 
implemented in a combined mwe-tdse solver to describe 
absorption and ultrafast pulse reshaping in a macroscopic 
medium. 

We used the TSDE-projection method to investigate 
the dependence of the xuv absorption cross section on 
the wavelength of the laser dressing at 10 12 laser 
intensity. We found complex multiphoton physics for 
800 nm light and shorter wavelengths. For longer, midin- 
frared wavelengths, however, we showed that the impact 
of the laser dressing in the Is — > 2p transition in helium 
can be described in terms of a A-type three-level model 
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previously used to describe eit for x-rays 7]. As in the 
earlier study, the transparency in helium is caused pre- 
dominantly by Autler-Townes splitting brought about by 
the strong one photon coupling induced by the dressing 
laser, in this case between the 2p and 2s states. We in- 
vestigated the macroscopic reshaping of an ultrafast XUV 
pulse resonant with this transition, for the case when the 
transparency is induced by an ir pulse which is substan- 
tially shorter than the XUV pulse. This means that the 
absorption properties of the helium atom change dynami- 
cally during its interaction with the xuv light. We found 
rich temporal reshaping dynamics in which the atoms 
both absorb and subsequently emit the xuv radiation 
in a process strongly influenced by Rabi oscillations be- 
tween the 2s and 2p states. This leads to an xuv pulse 
emerging from the macroscopic medium which has been 
shortened from 60 fs to 10 fs and whose peak intensity has 
increased by approximately a factor of two. The increase 
in the peak intensity, which results from the coherent 
population pumped into the 2p state before the dressing 
pulse arrives, is a truly dynamical effect which cannot 
be described in terms of a single absorption cross section 
only. 

Our results open up several possibilities for future 
research on ultrafast quantum optics. The control of 
xuv absorption by laser dressing of helium enables, 
for example, the possibility for post-generation ultrafast 
shaping of xuv pulses 0). And although we have con- 
fined ourselves in this work to the case of one xuv field 
with a dressing laser, it is a straightforward extension 
of the method to treat multiple xuv frequencies some of 
which some could be resonant with dressed atomic tran- 



sitions, and the complex interferences that would result 
from this 1, 3]. Also, we have used moderately strong 
IR fields that do not cause any excitation on their own, 
but the time-dependent treatment is not limited to these 
intensities. Using higher IR intensities will lead to gen- 
eration of harmonics in the non-linear medium. Har- 
monics with energies below and slightly above the ion- 
ization threshold, for which absorption dynamics plays 
the largest role, have recently attracted a lot of atten- 
tion, for instance as a source of vuv and xuv freq uency 
combs 51] or as a seed for free electron lasers [52]. The 
treatment of these processes necessitates using methods 
we have developed in this paper, since the atomic ab- 
sorption and emission properties will be changing on an 
ultrafast time scale. 
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